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ABSTRACT 

We examine the ability of photocvaporativc disk winds to explain the low- velocity components observed 
in the forbidden line spectra of low-mass T Tauri stars. Using the analytic model of Shu, Johnstone, & 
Hollenbach (1993) and HoUenbach et al. (1994) as a basis, we examine the characteristics of photoevapo- 
rative outflows with hydrodynamic simulations. General results from the simulations agree well with the 
analytic predictions, although some small differences are present. Most importantly, the flow of material 
from the disk surface develops at smaller radii than in the analytic approximations and the flow-velocity 
from the disk surface is only one-third the sound speed. A detailed presentation of observational conse- 
quences of the model is given, including predicted line widths, blue-shifts, and integrated luminosities of 
observable sulfur and nitrogen emission lines. We demonstrate that these predictions are in agreement 
with current observational data on the low-velocity forbidden line emission of ionized species from T 
Tauri stars. This is in contrast with magnetic wind models, which systematically under-predict these 
forbidden line luminosities. However, the present model cannot easily account for the luminosities of 
neutral oxygen lines in T Tauri stars. 

Subject headings: accretion, accretion disks — hydrodynamics — stars: formation 



1. INTRODUCTION 

It is now generally accepted that most young low-mass 
stars are born with circumstellar disks and that the evolu- 
tion of these disks holds a key to planet formation. Obser- 
vational evidence for circumstellar disks takes many forms, 
from direct imaging (e.g., Jayawardhana et al. 2002) and 
imaging of disk silhouettes (e.g., Bally, O'Dell, & Mc- 
Caughrean 2000) through indirect measurements such as 
infrared excess (e.g., Lada et al. 2000). These data provide 
an estimate for the lifetime of disks around low-mass stars, 
T'disk ^ 6 X 10^ yr (Haisch, Lada, & Lada 2001), requiring 
the presence of an efficient disk dispersal mechanism. 

Hollenbach, Yorke, & Johnstone (2000) considered a va- 
riety of disk dispersal mechanisms and concluded that vis- 
cous accretion of the disk onto the central star (e.g. Hart- 
mann et al. 1998), although dominating erosion of the 
inner disk, is incapable of removing the entire disk mass 
in the required time. Alternatively, photoevaporation of 
the disk, while providing an effective mechanism for disk 
removal at large radii > 10 AU, cannot remove the in- 
ner material [Shu, Johnstone, & Hollenbach 1993 (here- 
after SJH93); Hollenbach et al. 1994 (hereafter, HJLS94); 
Johnstone, Hollenbach, & Bally 1998]. Acting together 
viscous accretion and photoevaporation were predicted to 
efficiently remove the entire disk. Numerical calculations 
by Clarke, Gendrin, & Sotomayor (2001) and Matsuyama, 
Johnstone, & Hartmann (2003a) have shown that the com- 
bined effects of photoevaporation and viscous accretion are 
capable of dispersing disks on 10^ — 10^ yr timescales. Ad- 
ditionally, the formation of gaps within the gaseous disk 
during the dispersal era may place constraints on the evo- 
lution of planetary orbits (Matsuyama, Johnstone, & Mur- 



ray 2003b). 

The photoevaporation model for disk dispersal has been 
shown to fit the observational data well in the case of ex- 
ternal heating via nearby massive stars (Bally et al. 1998; 
Johnstone, Hollenbach, & Bally 1998; Storzer & Hollen- 
bach 1998). However, with the exception of a few cases 
(e.g., MWC349A), evidence for disk photoevaporation due 
to heating from the central star (SJH93, HJLS94) is largely 
circumstantial. Low- velocity ~ lOkms"^, blue-shifted 
forbidden lines observed in the spectra of T Tauri stars 
[Hartigan, Edwards, & Ghandour 1995 (hereafter HEG95)] 
may provide evidence for the 10"* K thermal disk winds 
expected in the photoevaporation model. The theoretical 
models of SJH93 and HJLS94 do not, however, include de- 
tailed hydrodynamic calculations and thus it has not been 
possible to test the model directly against these observa- 
tions. 

In this paper we take a first look at the hydrodynamic 
flow of thermally-driven disk winds powered by photoc- 
vaporativc heating in order to determine the observational 
consequences of the photoevaporation model. We begin in 
§2 with a brief description of how photoevaporation works, 
concentrating on photoevaporation due to heating from a 
star at the center of the disk. In §3 we describe how the hy- 
drodynamical simulations are performed and their results. 
The predicted observable properties of the simulations are 
presented in §4 and are compared with existing data. Fur- 
ther discussion and conclusions are presented in §5. 

2. PHOTOEVAPORATION OF CIRCUMSTELLAR DISKS 

Photoevaporation of disks by ionizing photons is a con- 
ceptually simple process to understand: high energy ra- 
diation (e.g., EUV, X-rays) from the central star ionizes 
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hydrogen at the disk surface producing a hot 10^ K ion- 
ized layer, similar to HII region conditions. Near the star 
this coronal zone above the disk is expected to be almost 
static and bound to the disk due to its location deep in 
the potential well of the star, however, at large radii from 
the star the thermal layer should be unbound, powering a 
thermally-driven disk wind. The critical length scale for 
the model is given in SJH93 and HJLS94 by 



To = 



(1) 



where M* is the mass of the star and Cs ~ lOkms"^ is the 
sound speed of the ionized gas. 

In order to determine the number density of hydrogen 
atoms at the base of the ionized layer SJH93 and HJLS94 
assume that for r < Vg the atmosphere is in hydrostatic 
equilibrium such that 

n(r, z) = no(r) exp [-zV2 H{rY] , (r < r^) (2) 

where z is the height above the disk and H{r) is the scale 
height; 

H{r) = r, (^^) ' , (r < r,). (3) 

The ionization of the material above the disk is not dom- 
inated by the direct passage of EUV photons from the 
central star but rather through the indirect EUV photons 
produced by recombination within the ionized layer. This 
effective scattering is strong because one in three recombi- 
nations goes directly to the ground state, producing a new 
Lyman continuum photon. The angular scale height, as 
seen from the central star, grows with distance from the 
disk center enabling the top of the ionized atmosphere to 
intercept direct EUV photons and redirect these photons 
through recombinations downwards (and upwards). Ana- 
Ij^ic analysis and numerical computation (HJLS94) allow 
for a determination of the base density power-law profile: 
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In the above a2 is the recombination coefficient to all states 
except the ground state and is the ionizing (EUV) pho- 
ton flux from the central star. The numerical analysis fixes 
the order unity constant to be C — 0.1. 

At radii larger than rg the ionized layer is assumed to 
flow at approximately the sound speed and thus the scale 
height takes the form 
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In this case the angular scale height is constant and the 
atmosphere directly above the disk no longer has a clear 
line of sight to the central star; therefore, the density at 
the base of the atmosphere decreases much quicker than 



within rg\ 
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The mass-loss rate due to the flow from both sides of the 
disk, as a function of disk radius, can then be computed 
directly: 



dr 



= A-Krno{r)csmii, {r > Vg), 



(8) 



where m// = 1.35 is the mean mass per hydrogen 
atom. 

Only two parameters are required to fully specify the 
photoevaporation model: and HJLS94 examined 
the effects of photoevaporation of the circumstellar disk 
around massive stars where the ionizing flux is set by the 
mass of the star $*(M,). SJH93, however, considered 
the effects of photoevaporation of the circumstellar disk 
around low-mass stars and concluded that only during the 
early T Tauri stage when the star disk accretion shock pro- 
duces enhanced EUV radiation would photoevaporation 
be important. Matsuyama et al. (2003a) self-consistently 
computed the strength of the star disk accretion shock 
during the dispersal of the disk and concluded that high 
ionization rates, > lO^'^s"^ are expected for the ini- 
tial ~ lO^yrs. Recently doubts have been raised about 
the possibility of ionizing radiation penetrating the thick 
protostellar jet, which lies between the star and the disk 
(Shang et al. 2002). Additionally, the accretion shock it- 
self may be thick to ionizing radiation (Alexander, Clarke, 
& Pringle 2004). In this paper we assume that a strong 
ionizing radiation field is present, either from the accretion 
shock, the chromospherically active star, or some other ra- 
diation generating process associated with accretion and 
the release of potential energy. 

A modified photoevaporation model for disk dispersal is 
applicable to situations in which the low-mass stars with 
disks are located near massive stars, such as in the Orion 
Trapezium. Johnstone et al. (1998) showed that in such 
cases the geometry for both the radiative transfer of the 
ionizing radiation and the hydrodynamic flow is much sim- 
pler and the model results provide a good match to the 
observations. However, it is clear that such external evap- 
oration is not applicable to every low-mass star. In the 
case of photoevaporation due to an isolated central star, 
detailed fits to the observational data have not been pos- 
sible due to the poor understanding of the flow charac- 
teristics. Yet there are clear strengths to such a model, 
inc;lu(iing the ability to disperse the disk on the observa- 
tionally determined timescale (SJH93; HoUenbach et al. 
2000; Matsuyama et al 2003a) and the ability to produce 
a thermally-driven wind with characteristic speed Cg sim- 
ilar to low-velocity, blue-shifted forbidden line value ob- 
served in T Tauri star spectra (HEG95). The photoevap- 
oration model also requires only a few free parameters and 
thus the observational characteristics should be well deter- 
mined. This aspect of the model is in direct contrast with 
the magnetic wind calculations often used to explain the 
low-velocity, blue-shifted forbidden line radiation (Cabrit, 
Ferreira, & Raga 1999; Garcia et al. 2001a, 2001b). 

3. HYDRODYNAMIC DISK WIND MODELS 

As discussed in §2, a detailed and complete picture of 
the thermal disk wind can only be attained by explic- 
itly solving for the hydrodynamic flow characteristics. For 
this purpose, we make use of the publicly available ZEUS- 
2D hydro dynamical code (Stone & Norman 1992a, 1992b; 



Photoevaporation of Circumstellar Disks 



3 



Stone, Mihalas, & Norman 1992). This code has been used 
in previous studies of disk outflows (e.g., Lee et al. 2001) 
and is a good choice because of its great flexibihty in terms 
of geometry and grid spacing. 

The effects of magnetic fields, radiation, radiative trans- 
fer, dust, and self-gravity of the gas are neglected in these 
models. We also neglect the effects of a stellar wind or an 
inner jet on the dynamics of the flow. A discussion of how 
this omission could potentially influence our results and 
conclusions is provided in §5. 

3.1. A spherical test: the Parker wind problem 

Before moving on to the problem of winds from pho- 
toevaporating circumstellar disks, it is instructive to con- 
sider a simpler problem to test the method. The Parker 
solar wind problem (Parker 1963) is ideal for this pur- 
pose because it is a trans-sonic outflow problem with a 
known solution. The Parker wind problem considers the 
outflow properties of a thermally-driven spherically sym- 
metric wind from a gravitating point-mass (such as the 
Sun). Assuming the wind is isothermal and barotropic 
(i.e., the pressure exerted on the gas is a function of den- 
sity only), then there is a well-deflned analytic solution 
with the wind passing through a sonic point (It^l/cs = 1) 
at the radius r/vg = 1/2. 

To set-up the problem in ZEUS-2D, a simple spherical 
coordinate grid is utilized. The radial coordinate, which 
is chosen to span r/vg = 0.1 — 1.0, is divided linearly into 
100 cells, while the polar coordinate, for which the range 
^ = 0° — 90° is used, is divided linearly into 60 cells. The 
flow is assumed to be azimuthally symmetric in ZEUS-2D. 
The boundary conditions of the simulation are all made to 
be reflecting except the outer radial boundary, which is set 
as an outflow boundary, and the inner boundary, which is 
continually updated such that the density remains con- 
stant along the inner edge (representing the source of the 
wind, i.e., near the surface of the star). The density at this 
inner edge was normalized to 1 while the remainder of the 
cells in the grid were set to an initial density of 10~^. The 
initial velocity within each cell is set to zero and the sim- 
ulation is allowed to run until it reaches steady-state. In 
practice, the velocity of the gas as a function of radius and 
time is tracked and we verify that the simulation converges 
rapidly to a steady-state solution identical (to within a few 
percent accuracy) to the analytic solution, with the sonic 
point occurring at r/vg = 1/2. Since the simulation setup 
just presented is very similar to the required setup for the 
disk wind models discussed below, the Parker wind anal- 
ysis provides a good test of their reliability. 

3.2. A disk wind test: power-law models 

Below wc explore winds originating from simple power- 
law disk models. Before proceeding, we emphasize here 
that these models are of academic interest only. Our pur- 
pose is to examine and imderstand the character of their 
solutions. This, in turn, provides a further test of our 
methods and also gives a basis for understanding the more 
realistic hybrid case considered below. 

3.2.1. Setup 

We now consider winds from photoevaporating disks. 
First, a description is given of the simulation setup in 



ZEUS-2D and how the models are initialized. 

At large radii the flow is expected to be dominated by 
pressure gradients in the diverging flow and should asymp- 
totically become radial. Therefore, a spherical polar co- 
ordinate system is selected for these simulations. Since it 
is not known a priori where the flow will become sonic, 
the radial coordinate spans a large range, r/rg = 0.1 — 20 
and is divided logarithmically into 200 cells. The polar 
coordinate is chosen to span the range 9 = 0° — 90° (thus, 
wc arc simulating just one quadrant of the atmosphere) 
and is divided linearly into 100 cells. By using this setup, 
we avoid excessive computation times but still have high 
enough resolution near the base, where the pressure gra- 
dients are largest, to accurately capture the dynamics of 
the flow. 

Similar to the Parker wind model described above, all of 
the boundaries are made to be reflecting except the outer 
radial boundary, which is an outflow boundary, and the 
disk base, which is updated after each time step to the 
original density distribution (representing replenishment 
from the photoevaporation of the neutral disk). The disk 
directly below the ionized atmosphere is assumed to be 
an infinitely flat, large (i.e., occupies the entire base of 
the grid) reservoir of neutral material. The base density 
in the ionized atmosphere is initially set-up as a simple 
power-law distribution with 

Mr) = "9 (~) ' 

where a = 3/2 or 5/2 is the power-law index, spanning 
the range of model conditions (i.e. Eqns. 4 and 7). 

The ionized disk surface is laid down on the equatorial 
axis {6 = 90°) and the density of the disk (and the above 
atmosphere) is normalized such that Ug = 1. The disk is 
assigned no velocity except for a Keplerian angular veloc- 
ity. As done for the Parker wind problem, the remaining 
cells, which comprise the 'atmosphere', are initially filled 
with a low-density (10~^ rig) gas and are assigned no ve- 
locity. The results of the simulation arc not at all sensitive 
to the exact value of this density. The simulations are al- 
lowed to run until they reach steady-state solutions. 

3.2.2. Disk wind test results 

Plotted in Figures la and 2a are steady-state stream- 
lines and velocity contours for the two test cases a = 3/2 
and a — 5/2. The streamlines (solid lines) are lines of 
constant angular momentum, j^, (which is a conserved 
quantity) and demonstrate how the ionized gas fiows from 
the disk. They are spaced in Figs. 1 and 2 according to 
the total disk wind mass-loss rate, Mdw,tot (see the caption 
of Fig. 1). A comparison of the two figures clearly shows 
that a much larger fraction of the total mass loss comes 
from the central regions of the a = 5/2 case than of the 
a = 3/2 case. Furthermore, the streamlines in the a = 5/2 
case arc approximately radial, whereas the streamlines in 
the a = 3/2 case contain a more significant vertical com- 
ponent. 

The differences in the spacings and directions of the 
streamlines of the two cases can be understood physically. 
Since the a = 5/2 case has a steeper base density profile 
than the a = 3/2 larger fraction of the ionized gas 

mass is located within the central region of that model. 
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Fig. 1. Steady-state results for the a = 3/2 case. Left: Streamline 
and velocity contours. The solid lines are streamlines and indicate 
how the ionized gas flows off the disk. They are spaced evenly ac- 
cording to the total mass-loss rate. Namely, the streamlines enclose 
20%, 40%, 60%, 80%, and 95% (not shown) of Md-w,tot- The dashed 
lines indicate where vtot equals 1, 2, and 3 times the sound speed 
(with the velocity of the flow increasing with distance from the disk) . 
Right: Density contours. Contours are evenly separated in log-space 
and range from logn/ug = —2.5 to 0. 
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Fig. 2. Steady-state results for the a = 5/2 case. The panels and 
lines have the same meaning as in Fig. 1. 

Indeed ~ 90% the total mass loss comes from inside rg for 
the a = 5/2 case. In addition, because the a = 5/2 case 
has a steeper density profile, this necessarily implies that it 
has a steeper pressure gradient as well. The steeper pres- 
sure gradient forces the streamlines to bend more quickly 
with respect to the vertical direction, producing a more 
spherical outflow. 

The dashed lines in Figs, la and 2a represent total ve- 
locity {vtot = yvf-\-vl) contours. The lines demonstrate 
where Vtot is equal to Cg, 2 Cg, and 3 c^. It is clear that 
the gas flows off the disk in the a = 5/2 case much faster 
than in the a = 3/2 case, directly related to the pressure 
gradients acting on the system. 

Intimately linked with the velocity and streamline prop- 
erties is the atmosphere density distribution. Plotted in 
Figs lb and 2b are density contours for the two cases. The 
contours are evenly spaced in log -space and range from 
logn/ug = -2.5 to 0. The a = 3/2 case shows horizontally 
elongated contours while the a = 5/2 case shows nearly 
spherical contours. Again, the difference arises because 
the a = 3/2 case has an appreciable mass- loss rate even 
at large radii, while a = 5/2 case has almost no mass loss 
at large radii. 

The small irregularities visible at large radii (near the 
base) in the a = 5/2 case are transient phenomena (not 
steady-state) and are due to the fact that the density of the 
base at these radii is comparable to that of the outflowing 
atmosphere and, therefore, the flow is easily disrupted. 
The low densities and mass-loss rates at these radii, 
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Fig. 3. Steady-state results for the truncated a = 3/2 case. The 
ionized disk has been truncated at r/vg = 5. The panels and lines 
have the same meaning as in Fig. 1. 

however, ensure that these irregularities have little or no 
effect on the rest of the (steady-state) atmosphere or on 
the predicted observable properties of these models. 

3.2.3. A truncated disk 

Figure 1 illustrates that a significant fraction of the to- 
tal mass- loss rate for the a = 3/2 case originates at large 
radii. This gives rise to the elongated density contours 
shown in Fig. lb. Furthermore, because of the shallow 
pressure profile of the ionized disk wind, the streamlines 
show large vertical components (i.e., perpendicular to the 
disk). We test these ideas quantitatively by running an- 
other a = 3/2 simulation but this time with the disk trun- 
cated at r/vg = 5. If the reasoning above is correct, we 
should expect the density contours for the truncated case 
to be more spherically symmetric while streamlines origi- 
nating near the truncation radius should show greatly re- 
duced vertical components (as there is no force at large 
radii pushing the gas vertically). 

Plotted in Figure 3 are the streamlines and velocity and 
density contours for the truncated a = 3/2 case. As ex- 
pected, the streamlines (Fig. 3a) illustrate that the flow is 
becoming approximately radial at large radii. Of note is 
the streamline which encloses 95% of the total mass-loss 
rate (and originates quite close to the truncation radius), 
which is essentially parallel to the ionized disk. Moreover, 
the density contours (Fig. 3b) also show a more spherical 
outflow than for the original a = 3/2 case plotted in Fig. 
lb. At large radii, the flow tends to spherical symmetry. 

3.2.4. Comparison to the Parker wind 

In the a = 3/2 case the disk mass loss is dominated 
by loss from the disk at large distances from the central 
star, producing a non-spherically symmetric outflow even 
at large radii. Alternatively, in the a = 5/2 case the disk 
mass loss is dominated by the central region and the far 
field limit tends to spherical symmetry. Comparison of the 
evaporating disk results with the analytic Parker wind so- 
lution provides an ideal framework for understanding the 
hydrodynamic flow. 

In Figure 4, the velocity and density profiles of the two 
cases are compared to the analytic solution to the Parker 
wind problem (assuming isothermal, barotropic gas). The 
simulation profiles are derived at a fixed polar coordinate 
9 = 45° but could have just as easily been chosen from 
any other slice through the atmosphere (except near the 
extremes 6* « 0° or 90°). 

The results plotted in Fig. 4a clearly indicate that the 
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Fig. 4. Velocity and density profiles at 9 = 45°. Left: Velocity 
profiles. Right: Density profiles. The short dashed line represents 
the « = 3/2 case, the long dashed line represents the a = 5/2 case, 
and the dot dashed line represents the PDW model (see §3.3). Also 
shown (solid line) is the prediction of the Parker wind model. The 
density profile of the Parker wind model has been arbitrarily nor- 
malized to match the other cases near r/rg ~ 2. For comparison, 
the dotted line in the right hand panel shows n oc r~^. 

shape of the velocity profile for the a = 5/2 case closely 
matches that of the Parker wind solution, including the 
normalization. The a = 3/2 case, however, has a much 
different shape and normalization than the Parker wind 
profile. This simply reflects the fact the a = 3/2 has a sig- 
nificant mass-loss rate even at large radii and the pressure 
forces from this gas are still modifying the dynamics of the 
flow (i.e., it is not yet dominated by spherical divergence). 

Fig. 4b illustrates that the flow from the a = 5/2 case 
is thermally-driven at large radii. For r > Vg, the density 
scales almost exactly as r~^, as expected for a thermally- 
driven flow (see the Parker wind line, for example). The 
a = 3/2 case, however, shows an over-dense flow at large 
radii with respect to the other cases. Again, this may 
be attributed to the non-negligible mass-loss rate at large 
radii in this case. 

3.2.5. Comparison to the analytic model 

The simulations behave as expected, at least at large 
radii, and we are now in a position to compare the steady- 
state properties with the model of SJH93 and HJLS94 
(from which we derive our physically-motivated initial con- 
ditions). As discussed in §2, hydrodynamic simulations, 
such as those performed here, are required for a proper 
treatment of the flow dynamics. Thus, the model of SJH93 
and HJLS94 incorporated a number of simplifying assump- 
tions about the form of the atmosphere, especially for the 
scale height H{r). In this section we test whether these 
assumptions are valid or not. 

In the SJII93 and IIJLS94 model, the density distribu- 
tion within rg, their "no flow" region, is determined en- 
tirely by the base density distribution, no(r), and by the 
scale height of the ionized atmosphere, H{r). Since we 
assume a base density distribution that is nearly identical 
to theirs, any differences between the models will manifest 
themselves in differences between the scale height distri- 
butions. Plotted in Figure 5 are the predicted scale height 
profiles for the a = 3/2 and a = 5/2 cases. Also shown 
is the analytic result of SJII93 and IIJLS94. Reassuringly, 
each of the cases closely follow the analytic result out to 
r/rg 0.4, matching both the shape and normalization 
of the analytic result. They each break away from the 
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Fig. 5. The scale height of the ionized atmosphere as a function 
of disk radius. The solid line is the analytic result of SJH93 and 
HJLS94. The short dashed line represents the a = 3/2 case, the 
long dashed line represents the a = case, and the dot dashed 
hne is for the PDW model (see §3.3). 

predicted relation at larger radii. This is not unexpected as 
the identification of Vg as the launching point for the out- 
flow is only accurate up to an order unity coefficient. For 
example, the Parker wind 'launches' at 0.5 rg due to pres- 
sure gradients in the divergent flow, and the disk wind can 
produce an even stronger pressure gradient. The fact that 
the simulations match the analytic result out to r/rg ~ 0.4 
implies that simulations and analytic model are neither 
hugely discrepant, nor are they strongly influenced by the 
base density distribution [at least in determining H{r)]. 
Thus, we are confident that a hydrodynamic model with 
a shallow base density slope near the star and a steeper 
density slope at large radii (such as that considered im- 
mediately below) will produce a reasonable model of the 
photoevaporating disks. 

3.3. Photoevaporative disk wind (PDW) model 

While the a = 3/2 and a = 5/2 cases were used to 
represent the range of base conditions, neither provides 
a reasonable fit to the full range of densities at the base 
of the photoevaporating disk models developed by SJH93 
and HJLS94. A better solution is to use a hybrid power- 
law model, which has a base density profile that mimics 
the a = 3/2 case at small radii and the a = 5/2 case at 
large radii. In fact, the radiative transfer calculations of 
HJLS94 demand such a model. We choose the following 
fitting function in order to quickly yet smoothly make the 
transition: 



noir) 



[{r/rgy/' 



(r/rg)25/2] 



1/5 



(10) 



Alternative fitting functions were tested and no signifi- 
cant changes in the results below were found. This hybrid 
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Fig. 6. Steady-state results for the photoevaporative disk wind 
model. The panels and lines have the same meaning as in Fig. 1. 

model, which, henceforth, will be referred to as the photo- 
evaporative disk wind (PDW) model, is nm in ZEUS-2D 
under the same conditions as the a = 3/2 and a = 5/2 
cases, as described in §3.2.1. 

Plotted in Figure 6 arc streamlines, velocity, and den- 
sity contours for the PDW model. A quick comparison 
of Fig. 6 with Figs. 1 and 2 reveals that, as expected, 
the PDW model is intermediate between the a = 3/2 and 
a = 5/2 cases. However, it more closely resembles the 
a = 5/2 case. The reason for this is that most of the total 
mass- loss rate in the PDW model is confined to r/rg < 2, 
which is very similar to the a = 5/2 case. In addition, 
both models appear to be approaching spherical symme- 
try at large radii (the density contours are nearly spherical 
while the streamlines are approximately radial). This is 
an expected result. Figs. 4 and 5 compare the results for 
the PDW model flow against the other two cases and the 
Parker wind solution. 

A physically insightful comparison between the PDW 
model and the; analytic model of SJH93 and HJLS94 is 
through the integrated mass-loss rate profile. The inte- 
grated mass-loss rate along the disk, assuming azimuthal 
symmetry and for both disk hemispheres, is calculated via 

Md^(r-) =47r/ p{r')ve{r'y dr' . (11) 

JO.lrg 

The mass-loss rate is presented in dimensionless units in 
Fig. 7. 

There are two noticeable differences between the PDW 
model and the analytic result. First, the PDW model 
shows significant mass-loss at radii smaller than Vg, 
whereas the analytic model has no mass-loss at these radii 
(by definition). Also clearly evident is that the PDW 
model predicts a lower total integrated mass-loss rate than 
the analytic result (by approximately a factor of 2.7). This 
is due entirely to the velocity of the flow off the disk being 
much smaller than that assumed by SJH93 and HJLS94 
(they assumed the flow velocity is equal to the sound 
speed). 

The total mass-loss rate is determined by both the base 
velocity and density. Since the models used here flx the 
density along the base to be nearly identical to the an- 
alytic result of HJLS94, it is the change in the launch 
velocity which accounts for the difference in the mass-loss 
rate. Our simulations have outflow velocities in the range 
of ~ 0.3 — QAcg. As a check, if wc assume vg = Cg and 
integrate (11) from r/r^ = 1 to 20 (to mimic the analytic 
result), we find a total mass-loss rate that is within 10% 
of the analytic result. 

The velocity profile of the flow is determined primarily 
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Fig. 7. Integrated disk wind mass-loss rates as a function of disk 
radius. The dot dashed line represents the PDW model, while the 
solid line represents the analytic result of SJH93 and HJLS94. 

from the divergence of the streamlines and depends on den- 
sity contrast, not the absolute density. Provided that the 
density drops as steeply as a = 5/2 in the outer disk, the 
base velocity should be reasonably well determined by the 
hydrodynamic simulation. The normalization of the den- 
sity, however, depends on the detailed balance between in- 
cident ionizing radiation, recombinations within the flow, 
and the liberation of new ions at the base. HJLS94 solved 
this balance equation explicitly, assuming a constant flow, 
and found that near r^, where the flow develops, recombi- 
nations dominate or are comparable to the production of 
new ions. Thus the base density varies smoothly through 
Vg with the power-law change reflecting the diminishing 
importance of overhead redirected ionizing radiation. This 
may need to be modified somewhat for the present model 
due to the acceleration of the flow and resultant drop in the 
number density of ionized hydrogen. Thus, the mass-loss 
rate presented here should be thought of as a lower-limit. 
Farther from the central star the new ions liberated at the 
base of the disk will eventually dominate the detailed bal- 
ance, resulting in a steeper density profile and even less 
evaporation. The density determinations of HJLS94 are 
strengthened when a lower base velocity is adopted. 

What are the implications of our simulations for the 
study of photoevaporating disks? A primary goal of 
HJLS94 was to sec whether the photoevaporation of neu- 
tral disks could replenish the expanding ionized gas in 
ultra-compact HII regions, thus explaining the paradoxical 
long lifetimes of such systems. In this context, the lifetime 
of an HII region is inversely proportional to the disk wind 
mass-loss rate. Since the simulations predict mass-loss 
rates that are within a factor of 2 — 3 of the HJLS94 model, 
our results do not significantly modify the conclusions of 
these authors. Rather, they reaffirm the importance of the 
photoevaporation mechanism proposed by HJLS94. The 
same situation holds for the dispersal of disks around low- 
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mass stars (SJH93; Clarke et al. 2001; Matsuyama et al. 
2003a). We do not expect major changes in the conclusions 
of those papers when these more accurate dispersal values 
are used. Only the details of how the disk evolves near Vg, 
where the hydrodynamic simulations show the strongest 
disagreement with the analytic theory, should be effected. 
Thus, the importance of a gap forming within the disk 
around and preventing planetary migration may need 
to be recalculated (Matsuyama et al. 2003b). 

After the work of SJH93 and HJLS94 there has been a 
number of published observational studies on low- velocity 
outflows in young stellar systems. Most notable of these 
is the T Tauri survey of HEG95. Below, we consider the 
observational consequences of our simulations for photo- 
evaporating disks around young, low-mass T Tauri stars. 
We then compare our results with the observational data 
of HEG95. This is the focus of §4. We show that the 
simulations generally give reasonable fits to a number of 
observed trends. 

4. OBSERVATIONAL CONSEQUENCES AND COMPARISON 
TO HEG95 

The most common way to constrain the properties of 
low-velocity outflows from T Tauri stars is through obser- 
vations of forbidden emission lines (typically, optical lines) . 
Commonly measured lines include [SIIJA6731, [SII]A6716, 
[NII]A6583, [OI]A6300, and [OI]A5577. Here, we explore 
the predicted proflles and strengths of such lines from our 
simulations and how these depend on, for example, incli- 
nation angle of the system and what region of the system 
is being observed (in the event the system is resolved). 
We also compare the model results with the observational 
data of HEG95. 

4.1. Rescaling for T Tauri systems 

Before line profiles and luminosities can be calculated, 

the simulations must be scaled from dimensionless units 
into physical units. The three fundamental quantities that 
must be rescaled are the velocity of the gas, the density of 
the gas, and the size of the system. 

First, the velocity of the gas is rescaled by considering 
that for a gas with a solar abundance of elements and a 
temperature of 10^ K, the sound speed, Cg, is ~ 10 km 
s-i. 

The physical size of the system is calculated by rescaling 
rg. Assuming a sound speed of 10 km s~^, the value of rg 
is given by 



1.3 X 10' 



M, 



cm 



(12) 



In what follows, we focus on systems with O.lA'/0 < < 
2Mq, similar to those observed by HEG95. 

Finally, we re-scale the density of the simulations. As 
discussed in SJH93 and HJLS94, the value of Ug is deter- 
mined by the mass of the star (or, cquivalcntly, rg) and by 
the rate of ionizing photons, $*, that are emitted by the 
star 

-3/2 / ^ \ 1/2 

We assume rates of 10''° s~' < < 10^^ s~', which is 
compatible with what is expected from T Tauri stars that 



n„«4xloV^V"7^J^) cm-3 (13) 



are accreting gas on to their surfaces (e.g., Matsuyama et 
al. 2003a). In §5, we discuss the validity of this assump- 
tion. 

The forbidden line luminosity, L, at a velocity u is given 

by 

L{u) = ^= / e ~^e{n,T)n{e,r)—dV , (14) 



n 



where n is the velocity shift along the line-of-sight due to 

the motion of the gas. e is the line emissivity per ion. n^/n 
is abundance of the atomic species from which the line 
originates, and dV = r^ sin 0d9d<j)dr. We have assumed 
that each cell in the simulation emits a Gaussian line pro- 
file, which is re-normalized by (27rc2)~i/2 to give unity 
when integrated. 

The velocity shift, /x, due to a particular cell's motion is 
given by 

/U = [{v0 cos 9 + Vr sin 6) cos cp — sin ^] sin i (15) 

— {Vr COS — vg sin 9) cos i 
where i is the inclination angle of the system. The disk is 
'face-on' when i = 0° and 'edge-on' when i = 90°. 

The line emissivity e{n,T) is calculated with the soft- 
ware package fivelevel, which was kindly provided by G. 
Mcllema (private commimication) . This is a key step since 
the lines we will focus on have markedly different critical 
densities. 

Finally, we assume that our simulations have a typical 
solar abundance of elements. We adopt the solar abun- 
dance ratios of Anders & Grevesse (1989). We also as- 
sume that all of the atoms are in the ionization state of 
interest. For example, when calculating the 673lA line, 
we assume 100% of the sulfur atoms in the disk wind are 
singly-ionized. 

Line profiles are constructed by integrating cqn. (14) 
over r = 0.1 - 20rg, 6* = 0° - 90°, 'and ^ = 0° - 360° (i.e., 
one hemisphere only, see discussion below) over a range 
of velocities. Since the simulations are 2-dimensional, we 
assume azimuthal symmetry in these calculations. Total 
integrated forbidden line luminosities and surface bright- 
nesses are calculated by integrating eqn. (14) over all ve- 
locities. 

4.2. Global properties: ionized forbidden lines 

The majority of observations of T Tauri systems to date 
are unresolved (e.g. IIEG95). The nature of unresolved 
observations is such that they can be used to directly con- 
strain only the global properties of these systems, such as 
the total forbidden line luminosity and the mean spec- 
trum of the system. Increasingly, however, high resolution 
instruments, such as the Space Telescope Imaging Spec- 
trograph (STIS) on board the Hubble Space Telescope, 
are resolving nearby systems (e.g., Bacciotti et al. 2000, 
2002). In such cases, it is becoming possible to measure 
radial properties, such as forbidden line surface brightness 
profiles and spatially-resolved spectra. Therefore, we are 
interested in examining both the global and 'structural' 
observational properties of the present PDW model. We 
examine here (and in §4.3) the observational consequences 
in the event the system is unresolved, while in §4.4 we 
focus on the observational properties of a resolved system. 

Plotted in Figure 8 are the predicted line profiles for 
the Nil and SII lines as a function of inclination angle for 
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Fig. 8. Predicted line profiles as a function of inclination angle for 
the PDW model. The lines have all been arbitrajrily normalized so 

that their peaks are equal to 1. 

the PDW model. The hne profiles are calculated assum- 
ing = IMq and = 10^^ s~^ The intensity of the 
lines have been scaled to make the peaks equal to unity in 
order to bring out differences in the shifts/shapes of the 
predicted profiles. 

Focusing on the upper left-hand panel, wo see that all 
three lines have been shifted to the blue by approximately 
10 km s~^. A blue-shift is expected since the disk is face- 
on, the outflow is moving away from the disk, and only 
one hemisphere was simulated. Had we simulated both 
ionized hemispheres, but with no neutral disk in between, 
we would observe no shift whatsoever. The presence of a 
thick neutral disk, however, acts as an occulting surface 
(with effectively infinite optical depth). Since genuine ob- 
servations also show blue-shifts similar to those plotted in 
Fig. 8 (e.g., HEG95), this implies that thick neutral disks 
are present in these systems. 

The two sulfur lines are shifted more to the blue than 
the nitrogen line (by ~ 2 km s~^). The reason for this is 
that the sulfur lines have lower critical densities than the 
nitrogen line. Thus, they arc more sensitive to the outer, 
low-density regions of the atmosphere, where the gas is 
flowing at a higher velocity. The widths (FWHM) of the 
three lines plotted in the upper left-hand panel range from 
w 26 — 29 km s~^, with the sulfur lines being shghtly wider 
than the nitrogen line. This can also be attributed to the 
lower critical densities of the sulfur lines. 

Examining the other three panels in Fig. 8, it is clear 
that as the inclination angle increases the line profiles be- 
come slightly wider and less shifted to the blue (see Table 
1). As the disk approaches the edge-on orientation the 
Keplerian nature of the disk is more observable and the 
radial outflow component is no longer solely flowing away 
from the observer. At i = 90°, there are equal amounts 
of gas flowing towards and away from the observer, hence 
there is no shift. 
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Fig. 9. Comparison of [SII]A6731 line profiles. The PDW model 
profile is calculated assuming = 10''^ s~^, M« = IMq, and 
i = 0°. 

How do the predicted line profiles compare with ob- 
served profiles? HEG95 found that the forbidden line emis- 
sion from T Tauri stars could be divided into two compo- 
nents, 'low- velocity' and 'high- velocity'. The high- velocity 
component correlates well with central jet diagnostics and 
shows extent in resolved observations. The low-velocity 
component, however, does not correlate with the central 
jet diagnostics and appears unresolved in the images. The 
low-velocity component appears to be associated with a 
disk wind (Kwim & Tademaru 1988; 1995; Hirth, Mundt, 
& Solf 1997) similar to the winds examined in the present 
study. Thus, we compare our results with the observed 
'low-velocity' component only. 

Plotted in Figure 9 is a comparison of predicted and 
observed [SII]A6731 line profiles. The predicted line pro- 
file is calculated for a face-on PDW model and assumes 
IMq and = lO''^ s'^ We note that changing 
these parameters (within reason) has only a small effect 
on the resulting line width and blue-shift. The observed 
profiles were taken from Fig. 14 of HEG95. All of the lines 
(both theoretical and observed) have been scaled such that 
their peaks are equal to 1. Thus, the plot should be used 
to compare widths and blue-shifts only. A comparison of 
integrated luminosities is given below. 

It is clear that the predicted profile is not a perfect 
match to the data. The observed lines show a large amount 
of scatter in the blue-shift of their peaks. However, we 
find it encouraging that the predicted peak of our profile 
falls in between and quite close to the peaks of the ob- 
served profiles. Our model predicts comparable, although 
slightly narrower, line widths than seen in the profiles of 
DQ Tau and CW Tau (UY Aur is clearly much wider and 
possibly contains two components). Thus, we would argue 
that the thermally-driven PDW model provides a simple 
explanation for the kinematics of the observed flows. It 
would be interesting to understand the slight difference in 
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predicted and observed line widths, however. As indicated 
by Table 1, an inclined system will give rise to increased 
line widths. This probably cannot explain the difference 
entirely, however. Below, we argue that a central gap in 
the disk may be necessary in order to explain the observed 
luminosities of the Nil line at 6583A. In terms of flux, a 
central gap would obviously give more weighting to the 
high velocity gas at large radii. As explained in §4.4, this, 
in turn, will give rise to increased line widths and could 
possibly explain this slight systematic difference. Alter- 
natively, a small boost from a central 'X-wind' could also 
give rise to increased line widths. 

Figure 10 presents the integrated line luminosities for 
the sulfur and nitrogen lines as a function of and 
Typically, the model predicts luminosities of a few times 
1O~^L0 up to a few times 1Q~'^ Lq, with systems that 
have a higher central stellar mass being more luminous. 
The trend of higher luminosities with higher central stel- 
lar masses (for a fixed value of is easy to understand. 
If the emissivity per ion, e, scales as n (as would be the 
case if critical densities were unimportant), we should ex- 
pect the luminosity to be independent of the mass since 
L cc e n dV oc n see eqns. 12-13). 

However, the sulfur and nitrogen lines have critical den- 
sities that are comparable to the densities found in the 
outflow, thus e docs not quite scale as the density and, 
hence, L cx where /3 > 0. This is confirmed by com- 
paring the Nil relation with the SII relations. Note that 
the Nil line luminosity does not depend as strongly on M* . 
The reason for this is that Nil has a higher critical den- 
sity than does SII and, thus, a smaller fraction of outflow 
has densities equal to or exceeding the critical density. It 
also explains the apparent bends in the L vs. relations, 
since the critical densities become less relevant for systems 
with higher central stellar masses (eqn. 13). The above 
arguments hold true for fixed values of but it is worth 
noting that will probably have some mass dependence 
(e.g., IIJLS94). Thus, the two parameters are not entirely 
independent. 

A comparison with the observational results of HEG95 
is also presented in Figure 10. For the data, we use the 
equivalent widths of the observed SII, Nil, and 01 for- 
bidden lines (in Table 4 of HEG95) to compute the ob- 
served integrated luminosities. We note that these equiva- 
lent widths were derived from the low-velocity component 
only (—60 km s~^ < u < 60 km s~^), which is appropriate 
for this comparison since we have not included jets in our 
hydrodynamic simulations. 

The top panels of Fig. 10 illustrate that the photoe- 
vaporative disk wind model predicts SII luminosities that 
are in excellent agreement with the observational data of 
HEG95. This is remarkable given that the observed sys- 
tems span more than an order of magnitude in M* and 
more than two orders of magnitude in L and that there 
was essentially no fine-tuning of the model parameters. 

The agreement with the observed Nil luminosities 
(lower panel of Fig. 10), however, is not as good. In partic- 
ular, our model tends to be over-luminous with respect to 
the data. Can the difference between the model and data 
for the Nil line be reconciled without spoiling the excellent 
agreement between the observed and predicted SII lumi- 
nosities? As noted before, the Nil line has a higher critical 
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Fig. 10. Predicted total line luminosities for the disk wind models. 
The hashed swath represents the predictions of the I'DW model. The 
low luminosity bound is set by assuming = 10^*^ s~^ and i = 0°, 
while the high luminosity bound is set by assuming = lO**^ s~^ 
and i = 90°. Solid squares are the observational data of HEG95, 
which is discussed in §5. Open triangles pointing to low luminosities 
indicate upper limits. 

density than the SII lines. Thus, it is more sensitive to 
the central regions of the disk wind than are the SII lines. 
A potential way of reconciling both the Nil and SII lines 
simultaneously is if the central region of the ionized disk 
atmosphere is either less dense than assumed by our PDW 
model or entirely removed from the system. HJLS94 con- 
sidered the effect which strong central winds might have 
on the structure of the ionized disk atmosphere and con- 
cluded that although the inner disk atmosphere might be 
severely limited, the outer evaporating disk would not be 
severely affected. Thus, it is possible that the absence of a 
strong low-velocity Nil component in the observations is 
evidence for the lack of an inner ionized atmosphere. 

4.3. Global properties: neutral forbidden lines 

We have purposely left the discussion of 01 forbidden 
emission lines until now. This is because accounting for 
the observed line strength of [OI]A6300 and [OI]A5577 is 
difficult in the PDW model since most of the oxygen in 
the wind should be ionized. Indeed, the similarity in the 
ionization potential of both hydrogen and oxygen, and the 
ensuing charge exchange between the two species, should 
ensure that the fraction of neutral oxygen is similar to 
the fraction of neutral hydrogen (Osterbrock 1989). Thus, 
only at the location where the disk wind crosses the ion- 
ization front might one expect to find significant 01 emis- 
sion. Given that the critical densities for these lines are 
even higher than for nitrogen, the line emission is domi- 
nated by the dense gas residing in the interior of the disk 
and near the base of the disk, just the location where the 
oxygen may be significantly neutral. Therefore, it is not 
immediately clear that the present model is completely 
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incapable of explaining the observed 01 emission. 

Indeed, externally evaporated disks are observed to pro- 
duce strong 01 emission (Bally et al. 1998) which is spa- 
tially confined to both the surface of the disk and the ion- 
ization front, with about equal total intensity in each com- 
ponent. The disk 01 emission must, therefore, arise within 
the photodissociation region (PDR) of the flow, below the 
zone where the ionizing (EUV) photons penetrate. 

Storzer & Hollenbach (1998) provide an explanation for 
this strong 01 emission from externally evaporated disk 
surfaces. Theoretical calculations of the dissociation of 
oxygen from OH (van Dishoeck & Dalgarno 1983) reveal 
that a significant fraction of the oxygen is left in an excited 
state. Storzer & Hollenbach determine the dissociation 
rate of OH and calculate the resulting 01 line luminosity. 
The enhanced far ultraviolet radiation field, due to the 
proximity of a massive star, and the high number density 
within the PDR, produce a strong 01 signature with an 
[O1JA6300 surface brightness of 2.2 x 10~^ ergss"^ sr~-^. 
For typical disk sizes of 100 AU this results in a total line 
luminosity of ~ 10""' L©. These values provide excellent 
agreement with the observations of Bally et al. (1998). 
However, because the 01 layer is confined to the base at- 
mosphere, where the outflow velocities are smallest, the 
resulting lines are likely to have only modest blue-shifts. 

Storzer & Hollenbach (2000) discuss the modification of 
the externally evaporating disk model for the case of disks 
evaporated by their central stars and determine that the 
far ultraviolet (FUV) radiation field is not strong enough 
to provide significant 01 line luminosities. This result, 
however, does not take into account the possibility of an 
ionized atmosphere redirecting EUV photons toward the 
disk. Indeed, the photocvaporating disk model of SJH93 
and HJLS94 requires the existence of a dense PDR sand- 
wiched between the ionized wind and the molecular disk. 

We check this idea quantitatively for the PDW model. 
The EUV photons come from a 10,000 K source and most 
will have energies near the Lyman limit of 13.6 cV. The 
cross section for hydrogen atoms to absorb these photons is 
of order 5 x 10~^* cm^. Thus, the EUV photons penetrate 
a neutral column of « 2 x lO"'^ cm^. This means that the 
ionization front, the region where the gas is still 10,000 K 
but has a significant fraction of neutral oxygen and hydro- 
gen, has a column density of H atoms of « 2 x 10^'' cm^. 
We have determined the zone at the base of atmosphere 
in the PDW simulation which contains this column den- 
sity for a range of stellar masses and ionizing photon rates. 
Typically, we find that the predicted blue-shifts are of or- 
der 3—4 km s~^, which compares quite favorably to the ob- 
servations of HEG95 (they find an average offset of ~ 5 km 
s^^ for neutral oxygen lines), but the integrated luminosi- 
ties are only of order lO"*" — lO~^iy0, which is about two 
orders of magnitude lower than observed (HEG95). Thus, 
without some kind of modification, the present model can- 
not easily accommodate the strength of the observed 01 
emission from T Tauri stars. 



4.4. Radial properties 

A potentially powerful test of the model, which will be 
possible within the coming years, is through comparisons 
with observed structural properties. Plotted in Figure 11 
are predicted nitrogen and sulfur line profiles as a function 
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Fig. 11. Predicted line profiles as a function of radius for the face- 
on (i = 0°) PDW model. The lines have all been normalized so that 
the peak of the — 20rg line (i.e., the entire disk) equals 1. 

of projected radius for the face-on PDW model. The lines 
are calculated assuming = IMq and = 10^^ s~^. 
Each of the lines have been normalized such that the peak 
of the line profile for the entire system, i.e., for — 20rg, 
is unity. Thus, in addition to showing how the shapes of 
the line profiles change as a function of projected radius, 
the plots indicate how the various regions of the disk con- 
tribute to the mean spectrum of the system. 

Focusing on the shapes of the profiles in the three pan- 
els, it is clear that as one goes further from the projected 
center of the outflow the profiles become more skewed. In 
particular, at large projected radii there are 'wings' evi- 
dent on the blue sides of the profiles. The development of 
line wings has a physical explanation. At small projected 
radii from the center, most of the flux from the system orig- 
inates from the dense ionized disk. Because the velocity of 
the gas flowing off the ionized disk is fairly constant as a 
function of radius, a nearly perfect Gaussian line profile is 
the result. Since most of the flux comes from the central 
regions of the outflow, this also explains why the mean 
spectrum of the entire system is almost exactly Gaussian. 
However, as one moves to larger projected radii the role of 
the atmosphere above the disk, in terms of flux, becomes 
more important. The gas in the atmosphere is flowing 
faster than the gas near the ionized disk (e.g.. Fig. 6a), 
hence, the Gaussian profile is skewed toward the blue. The 
increased importance of the atmosphere at large projected 
radii also causes the peak of the line profile to shift fur- 
ther to the blue and the FWHM of the line to increase. 
We demonstrate this quantitatively by fitting Gaussians 
to the line profiles. The results are presented in Table 1. 
While the shift in the peak is modest (~ 2 km s~^) in 
going from the mean spectrum of the system to just fo- 
cusing on the outer regions of the system, the increase in 
the line width is actually fairly substantial (> 15%) and 
potentially observable. 
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Fig. 12. Predicted forbidden line surface brightness profiles for the 
face-on PDW model. 

It is also evident in Fig. 11 tliat a larger fraction of the 
total flux from the nitrogen line comes from the central 
regions of the outflow than from the two sulfur lines. The 
reason for this trend is that the nitrogen line is more sen- 
sitive to the high density regions of the atmosphere than 
are the sulfur lines. 

Lastly, plotted in Figure 12 are the prcxlicted surface 
brightness profiles for the sulfur and nitrogcni lines. These 
are also calculated assuming M^, = IM© and = 10"*^ 
s~^. The surface brightness profiles of the three lines all 
exhibit the same general features: a central core and a 
power-law drop off at large radii. The size of the core is 
dictated by the critical densities of the lines, with the Nil 
profile being the most concentrated because of the higher 
critical density of nitrogen. At large radii, the profiles drop 
off as roughly Sune oc r~^. This can be understood phys- 
ically since Sii„e oc e n r oc r (e oc at large 
radii, low densities). The slight differences in the slopes 
of the three profiles are attributable to differences in the 
critical densities. 

5. DISCUSSION AND CONCLUSIONS 

The evidence for self-produced photoevaporative disk 
winds around low-mass stars is largely circumstantial. The 

primary reason for this has been the lack of detailed cal- 
culations on the outflow properties of these winds to com- 
pare with observations. Photoevaporative disk wind mod- 
els have, however, been shown to flt the observational data 
of externally- heated low-mass stars in the vicinity of mas- 
sive stars (e.g., .lohnstone ct al. 1998; Bally et al. 1998; 
Storzer & Hollenbach 1998) because the disk wind geom- 
etry and flow characteristics are much simpler. 

In this paper we use hydrodynamic simulations to com- 
pute the properties of the photoevaporative disk wind 
model and, subsequently, the predicted observational sig- 
natures. We show that the hydrodynamic simulations con- 
verge to steady-state solutions and that a number of gen- 



eral properties (e.g., scale height distribution, integrated 
disk wind mass-loss rate) agree well with the earlier an- 
alytic results of SJH93 and HJLS94. The outflows be- 
come nearly spherically symmetric and match the ana- 
lytic solution to the Parker solar wind problem at large 
radii, as expected. We further demonstrate that the pho- 
toevaporative winds produce blue-shifted ionized forbid- 
den lines with typical offset velocities of ~ —10 km s~^, 
widths of ~ 30 km s~^, and integrated luminosities of 
10~^ — 1O~*L0. These values compare favorably with the 
observations of HEG95 (see Figs. 9 and 10). The strength 
of the observed 01 emission is diflacult to account for with 
the present model, however. 

The photoevaporative disk winds examined here are in- 
tuitive and physically-simplistic. They rely on the avail- 
ability of = 10^°~^^ ionizing photons per second to be 
emitted by the central star disk accretion shock. Theoreti- 
cal arguments support such high levels of ionizing emission 
(e.g., Matsuyama et al. 2003a). Observationally, it is ex- 
tremely difficult to directly pin down the value of for 
low-mass T Tauri stars. HEG95 measured the mass accre- 
tion rates onto the surface of the T Tauri stars, however, 
and these, in principle, can be used to estimate the value of 
$». Given the mass accretion rate and the stellar radius, 
we calculate the accretion shock luminosity of the gas as 
it falls onto the star's surface via 



(16) 



where Mace is the mass accretion rate and i?* is the radius 
of the T Tauri star (see also Matsuyama et al. 2003a). Us- 
ing the inferred values of Lace a-nd assuming black-body 
radiation with a characteristic temperature of 10,000 K 
(Johns-Krull et al. 2000; Gullbring et al. 2000), we in- 
deed estimate ~ ]^q4o-42 g-i rpj^jg jgjj^jg credence to 

the photoevaporative wind hypothesis. 

Utilizing these ultraviolet photons, however, may be 
difficult. As mentioned above, both the protostellar jet 
(Shang et al. 2002) and the accretion column (Alexander 
et al. 2004) are likely to be optically thick to ionizing ra- 
diation, protecting the disk. Nevertheless, energy released 
within the accretion disk, and in the fast-moving jet, may 
produce significant ionizing radiation. Alternatively, the 
disk may be heated by another source such as FUV pho- 
tons (Johnstone et al. 1998, Storzer & Hollenbach 1998). 

The calculations presented here reveal that heating the 
disk surface to 10** K will produce a thermal wind with fiow 
characteristics similar to those observed around T Tauri 
stars. Additionally, if the mass-loss rate is high enough, 
comparable to the evaporating disk model, the forbidden 
line profiles of SII and Nil can be explained. Curiously, if 
the lO^'K wind is not ionized by EUV photons the resul- 
tant 01 forbidden line radiation also would be comparable 
to that observed. Thus, the observations of HEG95 seem 
to require a partially ionized 10^ K wind." 

Alternative models have also been proposed to explain 
the observed low- velocity component. In fact, the pre- 
vailing view is that T Tauri winds are magnetically-driven 
(e.g., Anderson et al. 2003). This is largely because strong 
magnetic fields seem to be the only way to power the 
coUimated, high velocity jets that are responsible for the 
high-velocity component observed in T Tauri spectra. In 
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addition, magnetic wind models are able to reproduce a 
number of the observed features of the high- velocity com- 
ponent (e.g., Shang, Shu, & Glassgold 1998; Garcia et al. 
2001a, 2001b). It is less clear, however, if the magnetic 
wind models can account for the low- velocity component, 
which was the subject of investigation in §4. Below, we 
briefly discuss how the predicted low-velocity components 
of magnetic wind models compare with their observational 
counterparts. We argue that the photoevaporative disk 
wind model proposed in the present study matches the 
current suite of observations at least as well magnetic wind 
models. 

Generally speaking, magnetic wind models may be di- 
vided into two classes: 'X-winds' and 'disk winds'. The 
primary difference between these two classes is the phys- 
ical extent of the magnetic field and, therefore, the loca- 
tion of the outflow launching radius. In X-wind models, 
the magnetic field is contained entirely within the very in- 
ner regions of the disk immediately next to the star. Disk 
wind models, on the other hand, assume that there is a 
large scale magnetic field present which extends to several 
AU. Unfortunately, it is extremely difficult to work out the 
exact observational predictions of such models in a self- 
consistent manner (Garcia et al. 2001a, 2001b). This is 
one advantage the physically-simplistic photoevaporative 
disk wind model has over the magnetic wind models. 

While it is difficult to calculate the precise observational 
characteristics of magnetic wind models, there have been a 
number of basic predictions made which are testable with 
current observational data. In terms of the X-wind mod- 
els, the predicted launching radius for the disk outflow is 
extremely close to the star. However, recent high reso- 
lution observations suggest that the low-velocity compo- 
nent probably originates at 'large' disk radii (e.g., Bac- 
ciotti et al. 2002; Anderson et al. 2003), suggesting that 
some other mechanism is responsible for the observed low- 
velocity component. In fact, irrespective of the model, one 
should expect the outflow velocity to be of order the escape 
velocity at the region where the wind originated. This is 
a strong argument that the low- velocity component must 
originate at moderate radii ('^ 1 — 10 AU). 

Magnetic 'disk wind' models, on the other hand, are able 
to produce both high and low- velocity components (e.g., 
Cabrit et al. 1999; Garcia et al. 2001a, 2001b), with the 
low-velocity component being launched from reasonably 
large disk radii. Furthermore, such models predict that 
the high velocity component is spatially extended while 
the low-velocity component is confined within the central 
(projected) regions of the system, which is in good agree- 
ment with observations. However, the problem with 'disk 
wind' models is that the level of ionization at large disk 
radii (i.e., at the launching radius) is too low to allow for 



strong magnetic coupling (e.g., Gammie 1996). Hence, 
such models tend to under-predict the luminosity of the 
low- velocity component by a substantial margin (e.g., Gar- 
cia et al. 2001b). As we have shown, however, the PDW 
model presented here (which contains no magnetic field) 
can naturally produce an outflow that matches a number 
of observed trends seen in low-mass T Tauri stars, so long 
as the photons can penetrate from the central star out 
to large radii (or if an external mechanism for producing 
EUV photons is present). It is possible, however, that both 
photoevaporation and a large scale magnetic held could 
be contributing to the low- velocity component, especially 
since the higher level of ionization due to photoevapora- 
tion should increase the effectiveness of the magnetic disk 
wind. This two component model could possibly explain 
the relatively large amount of scatter in the forbidden line 
proflles and luminosities and why some stars apparently 
show multiple low- velocity components (e.g., UY Aur, see 
Fig. 9). 

Additionally, it is likely that an interface between the 
fast jet-like wind arising from the inner disk will modify 
these results. HJLS94 modeled the effect of a strong stel- 
lar wind on the ionized atmosphere and concluded that 
only the inner regions would be affected. This suggests 
that the results of the calculations in this paper would 
need to be modified only in the central densest regions, 
resulting in somewhat lower line fluxes. This would pre- 
dominantly affect transitions having high critical densities 
such as the nitrogen line. Moreover, Matsuyama et al. 
(2003a) found that gaps could form in the inner part of the 
disk, near Vg, due to the erosive power of photoevapora- 
tion coupled to the viscous evolution of the disk. Although 
the Matsuyama models should be recalculated, taking into 
account the results of this paper, these gaps would remove 
the inner densest part of the ionized atmosphere produced 
in this paper. An evolutionary sequence might be visible 
in the low-velocity spectra from T Tauri stars dependent 
upon the state of the underlying disk. Detaik;(l calcula- 
tions, taking into account each of these factors, should be 
undertaken. 
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Table 1 



Line properties 
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FWHM 


shift 


FWHM 
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FWHM 
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i (deg.) 


(km s-i) 


(km s-i) 


(km s-i) 


(km s-i) 


(km s-i) 


(km s-i) 





26.5 


-9.2 


28.2 


-10.9 


29.4 


-11.8 


30 


28.7 


-8.1 


29.8 


-9.6 


31.0 


-10.3 


60 


33.2 


-4.9 


33.2 


-5.6 


34.2 


-6.0 


90 


35.4 


0.0 


34.9 


0.0 


35.8 


0.0 





FWHM 


shift 


FWHM 
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FWHM 
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Region 


(km s^"'") 


(km s""'^) 


(km s^""") 


(km s^"'") 


(km s^"^) 


(km s^ 


- 20 rg 


26.5 


-9.2 


28.2 


-10.9 


29.4 


-11.8 


1 - 20 rg 


28.0 


-10.3 


29.2 


-11.4 


30.4 


-12.3 


3 - 20 Tg 


31.6 


-11.8 


31.5 


-11.8 


32.5 


-12.7 


5 - 20 rg 


32.8 


-12.6 


32.7 


-12.5 


33.1 


-13.0 



